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INTRODU CTION 


Modern high-work turbines can be compact, transonic, supersonic, counter-rotating, or use 
a dense drive gas. The vast majority of modern rocket turbine designs fall mto these cate- 
gories. These turbines usually have large temperature variations across a given stage, and are 
characterized by large amounts of flow unsteadiness. The flow unsteadiness can have a major 
impact on the turbine performance and durability. For example, the Space Transportation 
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Main Engine (STME) fuel turbine, a high work, transonic design, was found to have an un- 
steady inter-row shock which reduced efficiency by 2 points and increased dynamic loading 
by 24 percent. The Revolutionary Reusable Technology Turbopump (RRTT), which uses full 
flow oxygen for its drive gas, was found to shed vortices with such energy as to raise serious 
blade durability concerns. In both cases, the sources of the problems were uncovered (before 
turbopump. testing) with the application of validated, unsteady computational fluid dynamics 
(CFD) to the designs. In the case of the RRTT and the Alternate Turbopump Development 
(ATD) turbines, the unsteady CFD codes have been used not just to identify problems, but to 
guide designs which mitigate problems due to unsteadiness. Using unsteady flow analyses as 
a part of the design process has led to turbine designs with higher performance (which affects 
temperature and mass flow rate) and fewer dynamics problems. The works of Griffin et al. [1], 
Garcia et al. [2] and Griffin and Dorney [3] are examples of the application of unsteady CFD 

to rocket turbine designs. 

One of the many assumptions made during the design and analysis of supersonic turbine 
stages is that the values of the specific heats are constant. In some analyses the value is based 
on an average of the expected upstream and downstream temperatures. In stages where the 
temperature can vary by 300 to 500°K, however, the assumption of constant fluid properties 
may lead to erroneous performance and durability predictions. In this study the suitability 
of assuming constant specific heats has been investigated by performing three-dimensional 
unsteady Navier-Stokes simulations for a supersonic turbine stage. In the first simulation the 
fluid properties were held constant, while in the second simulation the specific heats were 
allowed to vary locally as a function of static temperature. The predicted results from both 
simulations have been compared to the results of a well-validated meanline analysis which was 
also modified for variable specific heats. The meanline analysis is composed of one-dimensional 
equations of motion and empirical correlations. 

NUMERICAL ALGORITHM 

The governing equations considered in this study are the time dependent, three-dimensional 
Reynolds-averaged Navier-Stokes equations. To extend the equations of motion to turbulent 
flows, an eddy viscosity formulation is used. The turbulent viscosity is calculated using a 
two-layer algebraic turbulence model. The numerical algorithm consists of a time-marching, 
implicit, finite-difference scheme. The procedure is third-order spatially accurate and second- 


2 


order temporally accurate. The inviscid fluxes are discretized using a third-order upwind- 
biased scheme, while the viscous fluxes are calculated using standard central differences. An 
approximate-factorization technique is used to compute the time rate changes in the con- 
served variables. The numerical analysis uses Message Passing Interface (MPI) to reduce the 
computation time for large-scale three-dimensional simulations. 

The Navier-Stokes analysis uses 0- and H-type zonal grids to discretize the flow field and 
facilitate relative motion of the rotating components. The O-grids are body-fitted to the 
surfaces of the airfoils and generated using an elliptic equation solution procedure. They are 
used to properly resolve the viscous flow in the blade passages and to easily apply the algebraic 
turbulence model. The algebraically-generated H-grids are used to discretize the remainder of 

the flow field. 

The computational analysis has been validated on several supersonic turbine geometries 
(e.g., Refs. [3] and [4]). 

GEOMETRY AND GRIDS 

The supersonic turbine configuration, typical of those proposed for a reusable launch vehicle, 
consists of 21 vane airfoils and 52 rotor blades. In the current effort a 2-vane/5-rotor blade 
count approximation has been made. To keep the pitch-to-chord ratio constant, the vane 
airfoils were scaled by factor of 21/20 and the rotor blades were scaled by a factor of 52/50. 
The rotor tip clearance was set at a typical value of approximately 2.5% of the rotor height. 
Two simulations have been performed for the turbine, one with constant specific heats (Case 
1) and one in which the specific heats vary as a function of temperature (Case 2). 

The grid dimensions (number of passagesxixj xk) for the turbine simulations are presented 
in Table 1. Previous simulations suggested that 38 spanwise planes are adequate to resolve 
the predominant flow features for this type of geometry [3]. Figure 1 illustrates the midspan 
section of the computational grids. The total number of grid points used to discretize the 
turbine was 1,363,022. The average value of y + , the non-dimensional distance of the first grid 
line above the surface was approximately 1.5 for the airfoil surfaces and 3.5 for the endwalls. 

The simulations were run on seven 250 MHz processors of an SGI Origin2000. Each sim- 
ulation was run for 12.5 global cycles (in excess of one complete rotor revolution) at 25,000 
iterations per cycle, where a global cycle is defined as the time it takes for the five rotor blades 
to pass the two vane airfoils. The value of 25,000 iterations per cycle was chosen to resolve 
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all the expected frequencies of interest. Each iteration required approximately 16 seconds 
computation time on seven processors. The time periodicity of the solutions was determined 
by plotting pressure traces at different points along the airfoil surfaces. 

NUMERICAL RESULTS 

The turbine under consideration has an inlet Mach number of M 0 = 0.12, an inlet static 
pressure of 15.2MPa, and an inlet static temperature of T 0 = 1232 K. The rotor rotates at 
Q = 31,300 RPM, the Reynolds number (based on the inlet conditions and the rotor axial 
chord) is approximately 5 x 10 5 and the ratio of the rotor exit static pressure to vane inlet 
total pressure is P 2 /P t0 = 0.1875. The operating fluid in the turbine is a gaseous hydrogen 
mixture. In Case 1 a constant specific heat ratio of 7 = 1-354 was assumed, while in Case 2 
the specific heat ratio was varied locally according to the following equations: 

_ (1) 

' (C p /R) - 1 

C p /R = 01 + a 2 * T + o 3 * T 2 + a 4 * T 3 + a 5 * T 4 (2) 

where the coefficients for hydrogen are [5]: 

cti = 0.395086911? + 01 

a 2 = 0.10207987# + 00 

a 3 = 0. 131244661? — 04 

a 4 = -0.76649284# - 07 
a 5 = 0.34503763 E - 10 

for temperatures between 300 and 1000 K, and 

0l = 0.36440206# + 02 

a 2 = 0.546148011? — 01 

a 3 = — 0.16091 151J57 — 04 
a 4 = 0.21478497# - 08 

a 5 = -0.10131180^-12 

for temperatures between 1000 and 5000 K. The addition of the subroutine to calculate the 
local values of the specific heats resulted in less than a 5% increase in CPU time. 

Tables 2 and 3 contain time mean of the mass-averaged relative frame flow quantities pre- 
dicted by the Navier-Stokes simulation and meanline analysis for the constant and variable 
flow property simulations, respectively. The values presented for the Navier-Stokes simula- 
tions were time- and mass-averaged at the computational inlet (located approximately one 
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axial chord upstream of the vane leading edge), midway between the vane trailing edge and 
rotor leading edge and at the computational exit (approximately 1.25 axial chord lengths 
downstream of the rotor trailing edge). 

In general, the results from the Navier-Stokes and meanline analyses display good agreement 
in the vane passage, but exhibit differences in the rotor passage (see Tables 2 and 3). The 
meanline analysis predicts more flow turning and higher rotor exit relative Mach numbers, 
while the Navier-Stokes analysis predicts higher losses (associated with the endwall and tip 
gap flows). The differences between the meanline and Navier-Stokes solutions for the rotor 
are not surprising because the meanline code relies on empirical loss models for which there 
are limited data for supersonic turbomachinery flows. The total-to-total efficiencies predicted 
by the two analyses show fair agreement, whereas the predicted total-to-static efficiencies and 
work are in excellent agreement. 

Comparing the constant and variable fluid property data, the meanline results indicate 
insignificant changes in both the total-to-total and total-to-static efficiencies, while the Navier- 
Stokes analysis predicts a small decrease in the efficiencies. Both analyses predict a slight 
decrease in the work extracted by the turbine for variable specific heats. 

Time-averaged pressure distributions for the vane and rotor at 10%, 50% and 90% span are 
shown in Figs. 2 and 3, respectively. The vane loadings in both simulations are similar on 
the pressure surface and upstream of the throat on the suction surface. Downstream of the 
throat the case with variable fluid properties displays slightly lower pressures (which agrees 
with the Mach numbers in Tables 2 and 3) and increased loading. The results from Case 2 
indicate slightly more loading and higher pressures on the rotor at 50% and 90% span, while 
the loadings are similar at 10% span although the pressures are significantly lower in Case 2. 
The differences at 10% span suggest the nature of the flow in the hub region may be altered 
by the variable specific heats. 

Figures 4 and 5 illustrate a comparison of unsteady pressure traces near the leading edge of 
the rotor at 10% and 50% span, respectively. While the pressure traces between 25% and 75% 
span were similar in both simulations, there were significant differences (in both magnitude 
and frequency content) at 10% and 90% span (i.e., in the endwall regions). Interrogation 
of the frequency spectrum associated with the unsteadiness at the 10% span location of the 
rotor (obtained via Fourier decompositions of the unsteady pressure traces) reveals that Case 2 
contains significantly more unsteadiness at the vane passing frequency (approximately 10,000 
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Hz), as well as increased unsteadiness below and above the vane passing frequency (see Fig. 6). 
In Fig. 6 the frequencies for Case 2 were offset by 500 Hz to facilitate comparisons with Case 
1. In an effort to determine if the differences in the unsteady pressures at 10% span are 
confined to the leading edge of the rotor, additional traces were compared on both the suction 
and pressure surfaces. The traces showed similar behavior to the one near the leading edge. 
Thus, the unsteady pressures appear to be different in regions normally associated with strong 

secondary flows. 

CONCLUSIONS 

A set of unsteady three-dimensional Navier-Stokes simulations has been used to investigate 
the effects of varying specific heats on the performance of a supersonic turbine stage. One 
simulation was performed assuming constant specific heats, while one was performed in which 
the specific heats varied locally as a function of the static temperature. The results of both 
simulations have been compared with the output from meanline simulations. 

The flow quantities predicted by the meanline and Navier-Stokes analyses (for both con- 
stant and variable flow properties) show good agreement at the exit of the vane, but exhibit 
significant differences at the rotor exit (which is to be expected because of the lack of data 
and correlations available for supersonic turbines). There is excellent agreement between the 
total-to-static efficiencies and work predicted by both analyses. 

The Navier-Stokes results indicate increased pressure unsteadiness in the endwall regions 
for variable specific heats. In addition, the losses increased and the efficiencies decreased 
in the variable fluid property simulation. The meanline analysis predicted little change in 
the efficiencies for variable fluid properties, but similar to the Navier-Stokes results indicates 
reduced work. Thus, the current results indicate that it is important to include the effects of 
variable specific heats. 
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Table 1. Computational grid dimensions. 


Table 2. Case 1 time-averaged relative reference frame flow quantities. 
Table 3. Case 2 time-averaged relative reference frame flow quantities. 
Figure 1. Midspan section of turbine grid. 

Figure 2. Time-averaged pressure - vane; Case 1, - - Case 2. 

Figure 3. Time-averaged pressure - rotor; — Case 1, - - Case 2. 

Figure 4. Pressure history - rotor - 10% span - near leading edge; Case 

Figure 5. Pressure history - rotor - 50% span - near leading edge; Case 

Figure 6. Pressure decomposition - rotor - 10% span - near leading edge, 
Case 2. 
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Grid Type 

Vane 

Rotor 

0 

2x121x31x38 

5x121x21x38 

H 

2x86x41x38 

5x82x21x38 

Tip 

- 

5x121x16x5 

Total Points 

553,052 

858,370 


Table 1: 
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Navier-Stokes 

Meanline 

Variable 

Vane 

Rotor 

Vane 

Rotor 

Mi n 


1.15 

09 

1.22 

M 0 ut 

1.42 


Kfi£fl 

1.30 



mm 

15.2 

3.87 

Pout (Mpa) 

4.20 

2.89 

3.87 

2.89 

Ht m (kJ/kg) 

1787 

1622 

1787 

1619 

Ht out (kJ/kg) 

1786 


1787 

1619 

Pt in (MPa) 

15.4 

9.49 

mm 


Pt out (MPa) 

13.7 

6.36 


SB 

An (deg) 

0.0 

63.7 

0.0 

BB 

Pout (deg) 

73.3 

-61.6 

73.8 

-70.5 

7 

1.354 

1.354 

1.354 

1.354 

Reaction 

- 

0.095 

- 

0.073 

W (KJ/kg) 

- 

284.3 

- 

285.0 

Vtt 

- 

0.608 

- 

0.635 

Vis 

- 

0.449 

— 

0.450 


Table 2: 
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Navier-Stokes 

Meanline 

Variable 

Vane 

Rotor 

Vane 

Rotor 

M in 


1.16 


1.21 

Mout 

1.42 

1.10 

1.51 

1.31 

Pin (MPa) 

15.2 


15.2 

3.90 

Pout (Mpa) 

4.20 




Ht m (kJ/kg) 

1787 

1629 


1512 

Ht out (kJ/kg) 

1786 

UsfcUH 



Pt in (MPa) 

15.4 

KiM 

15.4 


Pt out (MPa) 

13.4 

6.29 

13.9 


An (deg) 

0.0 

63.6 

0.0 

72.5 

Pout (deg) 

73.1 

-62.0 

74.2 

-71.1 

Tin 

1.363 

1.380 

1.357 

1.383 

T out 

1.380 

1.381 

1.383 

1.386 

Reaction 

- 

0.092 

- 


W (KJ/kg) 

- 


- 


Vtt 

- 

0.602 


0.634 

Vts 

- 

0.447 

- 

0.450 


Table 3: 
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